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We study the smectic A-C phase transition in biaxial disordered environments, e.g. fully 
anisotropic aerogel. We find that both the A and C phases belong to the universality class of 
the " XY Bragg glass " , and therefore have quasi-long-ranged translational smectic order. The 
phase transition itself belongs to a new universality class, which we study using an e = 7/2 — ci ex- 
pansion. We find a stable fixed point, which implies a continuous transition, the critical exponents 
of which we calculate. 
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I. INTRODUCTION 

The effect of quenehed disorder on condensed mat- 
ter systems has been a widely studied topic for many 
years (THl], both for practical reasons (since disorder is 
always present in real systems) and fundamental ones. 
One of the fascinating results to emerge from the study 
of this problem is that even arbitrarily weak disorder can 
drastically affect the fundamental properties of many sys- 
tems. Among other effects, such weak disorder can de- 
stroy many types of long ranged order (e.g., ferromag- 
netic order in systems with quenched random fields Q), 
and it can radically change the critical behavior of many 
phase transitions |3|. 

Such effects have been found in, e.g., 
superconductors [l], charge density waves d, 0], Joseph- 
son junction arrays superfluid helium in aerogel^, 
and ferromagnetic superconductors [l^ . 

Some of the most novel and dramatic effects of 
quenched disorder are found in liquid crystals confined 
in random porous media [TTI. [T2j . These intriguing sys- 
tems exhibit a variety of exotic "Bragg Glass" phases. 
They also undergo unique types of phase transitions [Tsj. 
one of whi ch, the Smectic A to Smectic C (hereafter, AC) 
transition[lJ, |T5|, is the subject of this paper. 

In the high temperature phase {A phase), the nematic 
director n, which points along the axis of alignment of 
the constituent long molecules that make up the smectic 
material, and the normal to the smectic layers are 
parallel. In the low temperature phase (C phase), n and 
N tilt away from each other. 

Here, we wish to study how this transition changes 
when the smectic is confined in a quenched, random, 
anisotropic environment. This can be achieved exper- 
imentally by confining the smectic within aerogel. By 
changing the preparation of the aerogel, it should be 
possible to study the smectic A-C transition in different 
anisotropic environments. 

In previous work [isj . the nature of this transition 
was studied for a liquid crystal confined in a uniaxially 
stretched aerogel. In this paper we consider the transi- 



tion in a biaxial disordered environment, in which the 
aerogel is first stretched along one axis, then compressed 
in an arbitrary direction perpendicular to the stretchflG}. 
As illustrated in Fig. [1] we refer the stretched direction 
as the z axis, the compressed direction as the h axis, and 
the direction perpendicular to both z and h as the s axis. 
Wc denote the entire plane perpendicular to z (i.e., the 
s — h plane) as _L. 

We find that, in this geometry, the AC transition is 
a continuous transition belonging to a new universal- 
ity class. Unusually, the upper critical dimension, be- 
low which the transition exhibits non-mean-field critical 
behavior, is due = 7/2. We have studied this phase tran- 
sition using an e = ^ — d expansion, and find that there 
exists a stable fixed point, which implies a second-order 
phase transition with universal critical exponents. Fur- 
thermore, we find anisotropic scaling at the transition. 

One of these universal exponents is the usual order 
parameter exponent /3, defined via the temperature de- 
pendence of the tilt angle 9 between the local molecular 
axis (i.e., the nematic director n), and the normal to the 
smectic layers N. As T — > T^^, this angle, which is 
the magnitude of the order parameter for this transition, 
obeys the universal scaling law: 



9 cx {Tag - T)^ 



with 



1 2 

= e 



(I.l) 



(1.2) 



In the high temperature phase {A phase) , both the ne- 
matic directors n and the layer normals N lie, on average 
along z. In the low temperature phase (C phase), n and 
iV tilt simultaneously away from z in opposite directions 
in the s ~ z plane. Near Tac, the angles 0n{T) between 
N and z and 9n(T) between h and z are not independent, 
but rather have the same power-law temperature depen- 
dence as the total angle 9{T) = 9]y{T) — 0„(T) between 
N and fi. The ratio of 9n{T) to 9n{T) is a non-universal 
(negative) constant. 
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This geometrical constraint forbids azimuthal rota- 
tions of n and N in the C phase. As a result, there are no 
additional "Goldstone modes" associated with the break- 
ing of azimuthal symmetry in this anisotropic situation, 
in contrast to the isotropic caseflj, [l3|- Indeed, we find 
that both the A and C phases belong to the univers ality 
class of the random field XY model (RFXY) [H, [H [l|| . 

The correlation length of the transition is defined as 
the length beyond which the two point real space corre- 
lations of the order parameter c become small. Here c 
is the projection of N along s; its magnitude is there- 
fore proportional to the tilt angle 9n- As T — > Tag, 
the correlation lengths £,s,h,z along s, h, and z become 
very different: they diverge as power laws with different 
exponents: 
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To first order in e, we find 
1 

2^ 



1 + ^6 + 0(6^), 



(1.3) 

(1.4) 
(1.5) 
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The critical exponents most closely analo gou s to the 
exponent r\ in conventional phase transitions jl9l [20| are 
those associated with "anomalous elasticity" right at the 
critical point. Anomalous elasticity means that the elas- 
tic moduli characterizing the smectic are not, as in con- 
ventional materials, constants, but rather become depen- 
dent on the length scale A or wavevector q on which they 
are measured. More specifically, the bend modulus K 
that determines the free energy cost of bending the smec- 
tic layers diverges as (f — )■ according to scaling laws 



K{q) = q7''"--fK{qh/q^'\qz/qi^); 



(1.7) 



while the compression modulus B characterizing the re- 
sistance of the smectic layers to being squeezed closer 
together vanishes according to 



B{q) = qrfB{q,Jqi\qJq^')- 



(1.8) 



The anisotropy exponents (^h and are defined as = 
2 — ?7if/2 and Q = 2 — {rjg + rifc)/2. Our e-expansion 
results for (s,z quoted below show that neither of these 
exponents equals one. This implies that the system ex- 
hibits anisotropic scaling. In fact, since Ch Cz, there 
is completely anisotropic scaling between all three direc- 
tions: soft (s), hard {h), and along the stretch axis (z). 

The strength of the quenched disorder - specifically, 
of the random torques and random compression exerted 
by the aerogel on the smectic layers - is characterized by 
the disorder invariance At^c where the subscripts t, c re- 
spectively denote "torque" and "compression" . This also 
becomes wavevector dependent right at the AC transi- 
tion, with an anisotropic scaling form: 



To leading order in e = | — d, the anomalous exponents 
rjK, 11 B, and r]t are given by 



Vk = 


— e-HC'(e ) « 0.3, 


(I.IO) 
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^e + O(e2)«0.22, 


(i.ii) 


Vt = 




(1.12) 


Vc = 


2- ^e + 0(e2) w 1.55. 
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(1.13) 



The anisotropy exponents are: 
{riB + tik) 



Cz 



2 
1.74 
2 
1.85 



2-i^. + 0(.2) 



^ = 2-^. + 0(.2) 
2 27 ^ ^ 



(1.14) 
(1.15) 



Finally, we find that the specific heat exponent a is 
given by 



a=-e + 0{e^). 



(1.16) 



Since e = 1/2 is quite small in the physical dimension 
d = 3, we expect that our predictions for a, /3, and r/B,K,t 
are quantitatively accurate. 

In addition, there are several exact scaling relations 
between these critical exponents: 



/3 
a 

Vc 
Vz 



^ [4d - 10 + (5 - d)77K - - 277t] , (1.17) 



2 - ^ [4d - 6 + (3 - rf)?7A- - 2?/t - r/s] , 



2 + Tjt - iiB - Vk, 

CzVs, 
ChVs- 



(1.18) 
(1.19) 
(1.20) 
(1.21) 



AtAq)^q7"''^ft.Aqhlqi\qzlqi'). 



(1.9) 



As in the AC transition in uniaxially disordered environ- 
ments (l3| , a does not obey the usual anisotropic hyper- 
scaling relation a = 2 — Vg — id — 2)vh — Vz] this is due 
to the strongly relevant disorder. 

The behavior of the correlation functions near the AC 
transition is more complicated than in conventional phase 
transitions; in particular, it cannot be summarized by a 
single correlation length exponent v. This is in part due 
to the fact, mentioned above, that the transition exhibits 
anisotropic scaling. It is also due to the fact that the 
correlations within the A and C phases themselves are 
quite non-trivial. 

Since both the A and C phases belong to the RFXY 
universality class, we can use the results of reference [l3l 
to obtain the x-ray scattering intensity (which is a probe 
of the smectic translational order) far away from the crit- 
ical point. The most surprising feature of this is that 
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neither phase shows true long-ranged translational or- 
der; rather, they exhibit qtiasi- long-ranged translational 
order, but with a universal exponent [U, The ex- 

perimental signature of this quasi-long-ranged order is 
a power-law, rather than delta function, singularity of 
scattering as the Bragg peak is approached: 



/„ (Sq) cx [i6qsf + ahiSqhf + a,{Sq,f 



-,(1.22) 



where n is the index of the Bragg peak, and Sq ~ q — Gn 
with Gn the reciprocal lattice vector. In the A phase 
Gn = with a the layer spacing constant. In the 

C phase G„ = ^^{cosOnz ± sin^ArJ), where "±" re- 
flects that the tilting of the layer normals breaks two- 
fold symmetry, and ("— ") is taken if the tilting is 
along s (— s). ah^z are non-universal constants of order 
1. The exponent ~3+^55n jjj^piicg ^jj^^ In{Sq) diverges 
as a power law only for the first two Bragg peaks . 

In contrast, right at the AC critical point, the system 
only has short-ranged translational order; as a result, the 
Bragg peaks are broad. 

The X-ray scattering pattern near, but not at, the crit- 
ical point interpolates between these two limits, as illus- 
trated in Fig. m Sufficiently far away from the center 
of the Bragg peak, the line shape of I{Sq) versus 6q is 
qualitatively Lorcnzian squared. Close to the center of 
the Bragg peak, the line shape diverges as the power law 
given by Eq. (jl.22[) . The crossover between these two 
distinct patterns - or, equivalently, the width of the spike 
on the top of the line shape, is isotropic in scaling and 
strongly temperature dependent, 

Sq"" (xexp{-A\t\-^) , (1.23) 

where A is non- universal and of order 1, and t is the 
reduced temperature t = ^Tac'^ ' ^ ^® universal, and we 
can bound it above and below: 

.25<n< .65 . (1.24) 

The width of the quasi sharp spike decreases extremely 
rapidly as T ^ Tkci and vanishes completely right at 
T = Tac- Defining a correlation length ^ for the AC 
transition via ^ = q~^ , we see that ^ diverges like a 
stretched exponential as the AC transition is approached; 
hence, the usual correlation length exponent v = oo, in 
this sense. 

Another detailed experimental probe of our predictions 
is light scattering. Consider first an experiment in which 
the magnitude q of the wavenumbcr q is varied for fixed, 
generic direction q of the scattering wavcvector q. By 
"generic direction g " , we mean a direction of q for which 
no component of g is ^ 1. We find that in such an exper- 
iment, in both the A and C phases, the light scattering 
intensity I{q) plotted as a function of of q shows three 
distinct scaling regimes: 



q 

constant 



9 » , 

g^(T)«g«e,-i (1.25) 





FIG. 1: The geometry of the anisotropic disordered environ- 
ment. 



I(5q) 



l(§q)-|Sq|- 



8q^ 
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FIG. 2: The X-ray scattering intensity as a function of 5q, 
the deviation from the center of the first Bragg peak, near 
the critical point. 



where q^{T) vanishes as T — >■ Tag according 

q^{q}^\T-TAcW (1-26) 
where the universal exponent </> is given by 



f]B - Vk 



27 



0(e2 



(1.27) 



The first equality in (|I.27p is exact, while the second 
follows from the known e-expansions p.l0[) . (jl.ll[) . and 
P.13|) . for riK, r]B, and rjc. 

This prediction for I{cj) as a function of q for generic 
directions g of g is plotted in Fig. [S] 

More information can be obtained by restricting the 
scattering to the subspace of g with g/j = 0. In this 
case, we find that the light scattering exhibits different 
behaviors in the six different regions of the qh-qz , qh = 
plane illustrated in Fig. |4l 

In this figure, the various boundaries between regions 
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arc given by: 



FIG. 3: Plot of /(g), as given by equation p.25[) . versus the 
magnitude g of g for a fixed "generic" direction g of g, where 
"generic direction" is defined in the text. 
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FIG. 4: The six regions of different forms for the light scatter- 
ing in the g^ = plane. The boundaries between the different 
regions are given by equations (|I. 28111. 33|l , while the light scat- 
tering within each region is given by equation p.34|) . 
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Ss 1 

^qs, 
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and the light scattering in each region obeys: 



(1.28) 
(1.29) 
(1.30) 

(1.31) 

(1.32) 
(1.33) 



,2 U 



I{q) CX <^ 



qs \ tr 



(e 



[qz^zY 

2 
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region 2 

region 3 
region 4 
region 5 
region 6 



r 



(1.34) 



The two new characteristic lengths and are given by: 



Cs 
Cz 



^z I ^/v 



if 



(1.35) 
(1.36) 
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with 



= 3 - -riK 



(1.37) 
(1.38) 



We also find tliat tlie transition temperature Tag de- 
pends on the anisotropy of the disorder. Compared to the 
bulk transition temperature it is lowered by the stretch 
of the aerogel but raised by the compression. Henee, 
there is no a priori way to predict whether the net effect 
of the anisotropy is to raise or lower the AC transition 
temperature. 

In this problem there are two ways continuing the 
model from 3 dimension to higher dimensions. The 
e = 7/2 — d expansion uses the hard analytical contin- 
uation, in which we take the /i-space to be d — 2 di- 
mensional, while both the z- and s-spaces are kept one 
dimensional. It is complemented by the e = 4 — d expan- 
sion, which uses the soft analytical continuation, where 
both h- and z-space are one dimensional, and s-space is 
d—2 dimensional. The e expansion also leads to a stable 
non-Gaussian fixed point, which again implies a second- 
order phase transition. However, the numerical estimate 
of the exponents are not in good agreement with the e 
expansion in d = 3. The technical details are given in 
Appendix [X] 

The remainder of this paper is organized as follows. 
In Sec. ini we derive our model for the AC transition 
in biaxial disordered environments. In Sec. IIIIl we de- 
rive the renormalization-group recursion relations for the 
transition. In Sees. IVland llVl we calculate the thermo- 
dynamic critical exponents, and wavevector dependences 
of the elastic constants and disorder variances, respec- 
tively. In Sec. I VIII we determine the experimentally ob- 
servable correlation functions and make predictions for 
X-ray and visible light scattering. In Sec. IVIIIl we show 
that the transition is stable against dislocations in the 
smectic layers. 



II. MODEL 

In keeping with earlier work on the AC transition in 
clean 14 1 and uniaxially disordered [l3| systems, we as- 
sume that the local layer displacement u{r) and the local 
ncmatic director n(r) are the only important fluctuating 
quantities in this problem j23j . The complete Hamilto- 
nian, like Gaul, can be divided into three parts. The first 
part is obtained by a slight modification of the Hamil- 
tonian for the smectic A phase in isotropic disordered 
environments 113. Il3l. KLTl . It reads 



Since in our problem the liquid crystal is in an anisotropic 
environment, the Hamiltonian Eq. (jll.l[) need not be in- 
variant under rigid rotations of the liquid crystal. There- 
fore, the values of the coefficients B, g and w are not 
required to satisfy the constraint B = g = w, which they 
do in the isotropic caseflJjU^l- The _L and z components 
of the random field /i(r) in Eq. pi.ip respectively incor- 
porate the random torques and random compressions ex- 
erted on the smectic by the aerogel. We take h{r) to be 
Gaussian, zero mean, and characterized by anisotropic 
short-ranged correlations: 



(II.2) 

where 5'^ = 1 for i = j, and both i and j taking values 
in one of the d—2 hard directions, and zero otherwise; 
dfj — 1 for i = j = s, and zero otherwise; and Sfj = 1 for 
i — j ~ z, and zero otherwise. The function V{u — (j){r)) 
is a periodic function with period a, the smectic layer 
spacing. The field 4'{r) is a random length uniformly dis- 
tributed between and a with short-ranged correlations. 
Thus, V{u — 4>{r)) represents the random pinning of the 
smectic layers. 

The second part of our Hamiltonian describes the cou- 
pling between layer normals TV and directors n. The sim- 
plest such coupling is given by [l^] 



H. = 



d'^1 



P\N 



(11.3) 



This part of the Hamiltonian is responsible for the AC 
transition. To see this, simply note that when P > 0, 
the lowest energy state has N \\ n; this corresponds to 
the A phase. In contrast, when P < 0, the energy can be 
lowered by making N and n point in different directions; 
this corresponds to the C phase. Thus, the point P = 
corresponds, in mean field theory, and in the absence of 
other couplings between N and h, to the critical point 
for the pure smectic AC transition. If the coupling P is 
assumed to be a monotonically decreasing function of the 
temperature T, the A to C transition will occur as the 
temperature is lowered through the temperature at which 
P = ( again, in mean field theory and the absence of 
other couplings between N and h). 

Since the layer normal N is related to the smectic layer 
displacement u by the geometric relation 



N : 



(11.4) 



wc can rewrite Eq. pi.3l) as a coupling between u and n; 
to lowest non-trivial order in u, this reads: 



Hi = I d^r 



^(Vi^^)^ + |(5,^.)2-|(a,«)|Vxup 



+ h{f)-Vu + Vpiu - 0(f)) 



(H.l) 



Ho 



d^r P\y±u + 6n±\ 



(H.5) 



Since the aerogel is stretched and compressed, the dis- 
order is no longer isotropic. Instead, on average it picks 
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out a preferred orientation for both the layer normals 
and the directors. Therefore, in general the Hamiltonian 
should also include 



Q\SnA 



+ (M' - M){dhuf + (Q' - Q){6nh)^] , (II.6) 

which is not rotation invariant and hence forbidden in 
the isotropic case. It is expected, that the stretching of 
the aerogel tends to hold both N and n along z, while 
the compression tends to hold both N and h within the 
zs plane. Therefore, the coefficients in Eq. (|II.6P obev 
the bounds 



M' > M >0 



Q' >Q>0 



(II.7) 



The values of M and Q reflect how much the aerogel is 
stretched; those of M' — M and Q' — Q likewise reflect 
the degree of compression. 

After a linear change of variables 



Sn' 



Sus + RdgU , 
Snii + R'dhU 



(II.. 



the sum of H2 and becomes 
■D{T) 



{d,uy + {P + Q){5n'^y 



+ l{d,uf + {P + Q'){Sn'^f 



(11.9) 



where R = P/{Q + P), R' = P/{Q' + P), D(T) = 
2M + 2QP/{Q + P), 7 = 2A/' + 2Q'P/{Q' + P). It 
is easy to verify that under the constraints (|II.7|) D{T) < 
^<P + Q<P+Q'. Therefore, if we start in the high 
temperature (A) phase and lower the temperature, D{T) 
is the first among the coefficients in Eq. (|II.9[) to become 
negative. As we will see in a moment, the critical point 
(ignoring fluctuations) is given by D{Tac) = 0, which 
leads to P{Tac) = -MQ/{Q + M). Therefore, stretch- 
ing the aerogel lowers the critical temperature. Stretch- 
ing aerogel suppresses the transition. Since Pac has no 
dependence on M' and Q', Tag is not affected by the 
compression, in this mean field approximation in which 
fluctuations are neglected. 

Now we look for the ground state of the system, av- 
eraged over many realizations of the disorder. Since 
the disorder pieces in Eq. (jll.l[) are completely random 
and isotropic, they have no averaged contribution to the 
ground state. Therefore, they can be neglected. Mini- 
mizing the sum of the Hamiltonians Eq. pi.l[) and Eq. 
pi.gp we find two distinct states separated by D{T) = 0. 
For D{T) > (i.e. T > Tag, ), the ground state is 
given by |V_lu| = \dn±\ = 0, which implies that both TV 
and n point along z. Thus the system is in the A phase. 
For D{T) < (i.e. T < Tag), the ground state shifts 
to dsu = ±y/-2D{T)/{w- g'^/B), Sn, = -Rd^u, dhu = 
Srih — 0. This implies that both N and n tilt away from 



z. However, the tilting is constrained since N and h are 
restricted within zs plane with z in between the two (this 
constraint arises from the fact that the last two terms in 
|II.9| will raise the energy of any state in which N and h 
do not lie in this plane above that of an otherwise iden- 
tical state in which they do). In this case the system is 
in the C phase. Near the critical point, the magnitudes 
of the two tilting angles 0jv and On are very small such 
that approximately, \9m{T)\ = \dsu\ and \9n{T)\ = \5ns\. 
Therefore, they obey the relation 



IN 



{T) P + Q' 



(11.10) 



Although Eq. III. 101 is derived through mean field the- 
ory, since we are simply minimizing the free energy while 
ignoring fluctuation effects, it still holds near the critical 
point. This is because the fluctuations of Sn'± are al- 
ways massive. For the same reason we can set 5n' ±^ = 
in (jll.9|) and use the sum of (jll.ip and (jll.9[) as our sim- 
plified model: 



H 



+ h{r)-Vu + Vp{u-(t>{r)) 



(11.11) 



where the anharmonic terms involving dhU are not in- 
cluded since they are irrelevant compared to ^{dhuY /2. 

To cope with the quenched disorder we employ the 
replica trick We assume that the free energy of the 
system for a specific realization of the disorder is the 
same as the one averaged over many realizations. To 
calculate the averaged free energy F = InZ, where Z 
is the partition function, we use the mathematical iden- 
tity In Z = lim„ 



^0 „ ■ When calculating Z", we can 

first compute the average over the two random fields h{r) 
and 4){r), whose statistics have been given earlier. Imple- 
menting this procedure gives an replicated Hamiltonian 
with the effect of the random fields transformed into cou- 
plings between replicated fields: 



a, 13 = 1 



At 
T 



Ac 1 
-7jT{dzUa){dzUj3) + —Ap{Ua - Up) 



(11.12) 



We will now analyze this Hamiltonian using the renor- 
malization group. 
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III. RENORMALIZATION GROUP 

Right at the critical point, since the tilt term {dsUa)^ 
vanishes, the Hamiltonian Eq. pi.l2[) becomes very simi- 
lar to that for a smectic A liquid crystal confined in com- 
pressed aerogel; such a system has been described [l7| as 
an "m = 1 Bragg Glass" . For the latter Ref. showed 
that under the renormalization group (RG) , the random 
pinning Ap is irrelevant, in the RG sense, in spatial di- 
mensions d < |- Repeating the (virtually identical ) cal- 
culation we find that within the critical region, Ap is also 
irrelevant in our problem. Furthermore, it is easy to show 
that the random compression Ac is also subdominant to 
At. Therefore, to construct the RG for the transition we 
can use the following truncated Hamiltonian: 



At 



5Z 

a, 0=1 



(III.l) 



where we have adopted the hard analytical continuation 
with a single stretch (z) axis, a single soft(s) axis, and 
d~ 2 hard (h) axes. 

First we calculate the harmonic propagator. In fourier 
space the quadratic part of the Hamiltonian can be writ- 
ten as 



H 



2 ^ 

9 



Aai3{q)Ua{q)up{-q} 



(111.2) 



where 

Aal3 

and 



(Bql+jql + Dq^+Kqt) S^^ 



At 2 



Ua{q) = 



(HI.3) 



(HI.4) 



The correlation function of this model can be shown to 
be 



{uo,{q)up{q')) ^TA^l{q)5i;^^^, 



(HI.5) 



where A^^ is the matrix inverse of A^p . For any n x n 
matrix M„b of the form 



MaB = aSaft + b. 



it is straightforward to verify by explicit matrix multipli- 
cation that 



M 



— SaB 

a 



a{a + bn) 



(111.7) 



Since Aq,/? is if the form (jHI.6p . with 
L>g2 + Kqj and 6 - ^ 
taking the limit n — > 



Bql 



111 



^q±, we obtain, from (jHI.Sp . upon 



{Uo.{q)up{-q)) = TGiq)S^p + AtqlG{q)^ (IIU 



with 



Bq^, + iqi + Dqj + Kqj 



(HI.9) 



The momentum shell RG procedure consists of trac- 
ing over the short wavelength Fourier modes of Ua{r) 
followed by a rescaling of the length. As usual in 
the momentum shell RG procedure, we initially restrict 
wavevectors to lie in a bounded Brillouin zone; here, we 
will take our Brillouin zone to be the infinite slab —A < 
qs < A, — oo < |(f;i| < oo, — oo < qz < oo, where A ^ 1/a 
is an ultra-violet cutoff, and a is the smectic layer spac- 
ing. The displacement field Ua{r) is separated into high 
and low wave vector parts Uaif) = W^if) + (r), where 
u„ (r) has support in the wave vector space Ae^"^^ < qs < 
K, —GO < qh < oo, —oo < qz < oo. We first integrate out 
u„ (r). This integration is done perturbatively in the an- 
harmonic couplings in IHI.ll as usual, this perturbation 
theory can be represented by Fcynmann graphs, with the 
order of perturbation theory reflected by the number of 
loops in the graphs we consider. After this perturbative 
step, we anisotropically rescale lengths, with = r^e^, 
rh = r-^e"''^, Tz = r'^e'^'^, so as to restore the UV cutoff 
back to A. This is then followed by rescaling the long 
wave length part of the field with u^{r) = u'^{r')e^^. Af- 
ter this procedure we obtain the following RG flow equa- 
tions to one-loop order, keeping only lowest order terms 
in D, since we are interested in the critical point where 
D, to leading order in e, vanishes: 



m 

dB 



dK 

m 

dg 



dw 



[{d - A)oJh +0Jz + 2x + l]j , 
{d - 2)ujh -ujz + ^x+l 



(HI.IO) 



32V2 



.91 



B 



(d - 2)ujh + ujz + 2x-i + 



8V2 



91 



(HI.ll) 
K, (HI.12) 



{d - 2)LUh + 3x - 1 + 



16V2^^ 32V2 



92 



(HI.13) 



{d - 2)L0h + + 4x - 3 - 



8^2 .92 



+ 



8\/f^^ 32V2 



27 



.92 



(HI.6) — 



dAt 
dl 
dD 



{d - 2)ujh + a;^ + 2x - 1 + 



[d - 2)L0h +uJz + 2x-l + 



91 



32V2. 



(HI.14) 
At, (HI.15) 



32v^ 



(51 - .92) 



D 



8^2 



K{g2 - 51) 



(HI.16) 
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gi and 52 are two dimcnsionless couplings defined by 
51 = Cd-i{g/B)^AtiB^'-''K''-y/^A'^'', (III.17) 



92 



= Cd-iiw/B)At{B-f^-'^K'^~y/^A^'^-'^. (III.18) 



where Cd = Sd/i^n)'^ with Sd the surface area of a d- 
dimensional sphere of radius one. The RG recursion re- 
lations pil.lll IIII.121 IIIL131 IIII.151 IIIL16|) imply that gi 
and 52 themselves flow according to the recursion rela- 
tions 



dgi 
M 

dg2 
M 



2e 



2e 



21 

64 
63 



91 



16^/2 
27 



64 v^^^ 32 v^'''^ 8V252 



.92 



.92 



91 1 



9l 



(III.19) 

52, 

(III.20) 



where e = 7/2 — d. These two flow equations are inde- 
pendent of the rescaling factors oj/j ^ and x since 51 2 are 
dimcnsionless. 

To discuss the critical behavior it is convenient, but 
not necessary, to make a special choice of loz^h such that 
7, i?, and K are fixed at their bare values. This choice 
is given by 



= 2- 

X 



1 



16V2'^' 
7 



b-2d 4(i - 9 



2 128%/2 
After this choice Eq. piI.16p becomes 



51- 



(III.21) 
(III.22) 
(III.23) 



dD 



5ffi _ 9g2 
32\/2 32V2. 



D 



8^2 



A'(g2 



For d < 7/2, there are three fixed points to Eqs. (IIII.191 
IIII.201 111124)) . All three lie on the "critical surface" [24} . 
which separates flows that asymptotically run to large, 
positive D from those that run to large, negative D. More 
physically, this surface is the phase boundary between the 
A and the C phase. The RG flow on this critical surface, 
projected on to the gi 2 parameter space, is illustrated in 
Fig. El 

This flow is topologically identical to that for the AC 
transition in a clean system in an external field in spatial 
dimensions (i = 3 — e, fore<Cl [13 ■ This is scarcely sur- 
prising, since our problem is simply a disordered version 
of that problem. 

The Gaussian fixed point G, given by g\ — 0, 52 = 
0, D* = is trivial, and unstable. Hence, it does not 
control the AC transition. 

The uniaxial fixed point U, given by g^ = 32 ~ 
0, is stable only within the pa- 
: 52, = 0. This fixed point was 



128\/2 



rameter subspace gi 



first found in the study of smectics in uniaxially com- 
pressed aerogel [iTj. and controls the phase reference (iTj 
calls the "m = 1 Bragg Glass" . In that system, the con- 
straint gi ~ g2, D = Q is enforced by the symmetry that 
the Hamiltonian is invariant under any rigid rotation of 
the liquid crystal about the compressed direction. This 
fixed point controls the power law anomalous elasticity 
of the m = 1 BG. 

In our problem, since the rotation invariancc is de- 
stroyed by the stretching of the aerogel, the conditions 
.91 = 92 and D = are not generally obeyed. Thus, we 
expect that the critical behavior is controlled by a differ- 
ent fixed point with g\ 7^ 1721 which we call the biaxial 
fixed point, and denote as B. In fact, the fixed point B 
must satisfy 5i < ^2 : as required by the stability of the 
model Eq. pi. lip . To see this we rewrite the model as 



H 



d^r 
D 



+2 



{dsu) 



2 , 1 ,c ..n2 1 



-{dhu)- + - u) 



2 ' ' 2 
+h{r)-Vu + Vp{u~(i){f)) 



B 



(III.25) 



Clearly, for the Hamiltonian to have a well defined mini- 
mum at It = the coefficient of (9sm)^ must be positive, 
which implies gi < (72- This stability condition should be 
preserved under the RG and, indeed, satisfied by fixed 
point B. To the first order in e, the fixed point B is given 
by 



9i 

92 

D* 



128\/2 
^7 
512\/2 
81 



e + 0{e% 



(IH.26) 
(IH.27) 
(IH.28) 



It is easy to show that this fixed point is stable in two 
directions (to leading order in e, the 51^2 plane), and un- 
stable in a single direction (to leading order in e, the 
D direction). This stability/instability structure is pre- 
cisely that associated in the standard theory of critical 
phenomenon [1^ with a fixed point controlling a phase 
transition. Therefore, the fixed point B controls the crit- 
ical behavior of the AC transition. 



IV. ANOMALOUS ELASTICITY WITHIN THE 
CRITICAL REGION 

It has long been known that both pure[iu| and 
disordered [l^ smectics exhibit anomalous elasticity, 
which is a novel phenomenon found, thus far, only in 
liquid crystalSjJuctuating membranes, and at phase 
transitions [1^ [g^]. In smectics, the anomalous elastic- 
ity is characterized by wavevector dependent compression 
modulus B and bend modulus K at large length scales. 
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by 




FIG. 5: Graphical RG flow of the dimensionless couphngs 
gi and g2, in the D = plane, for rf < |. The fixed point 

U was previously found in and controls 

m = 1 BG phase. Both the Gaussian fixed point G (0, 0) 
and the uniaxial fixed point U are unstable. The biaxial fixed 
point B (i2|^e, ^^j^e) is stable, and controls the AC tran- 
sition in biaxially disordered media. 



The origin of this exotic behavior is the non-linear cou- 
pling of fluctuations. For a smectic in an isotropic envi- 
ronment (either clean or disordered) , rotation invariance 
symmetry, which is spontaneously broken, prohibits cer- 
tain terms (e.g., | V_LMp) in the elastic model, and, hence, 
the system is soft in the _L directions. This softness leads 
to big fluctuations. As a result, the anharmonic terms in 
the model become nontrivial and the conventional har- 
monic theory breaks down. In the presence of quenched 
disorder, since the fluctuations are bigger, the anhar- 
monic effect becomes even more important. Thus, in 
contrast to the marginally weak anomalous elasticity in 
pure smectics with logarithmic wavevector dependence, 
disordered smectics show a strong one with power law 
wavevector dependence. 

In our problem, since rotation invariance is broken, in 
general all _L directions are hard. However, near the crit- 
ical point (i.e., at length scales shorter than the correla- 
tion lengths), when {dsUaf' is subdominant to {plua) 
and hence negligible, the softness in the s direction is 
recovered. Indeed, in this case the elastic theory of our 
system becomes very similar to that for the m = 1 BG 



phase of a smectic in uniaxially compressed aerogel [17 . 
Therefore, we expect similar anomalous behavior. How- 
ever, the anomalous exponents will be different, since 
the system flows under the RG to a different stable fixed 
point, as illustrated in Fig. [S] 



To get an estimate of how important the anharmonic 
effect is, we calculate the graphical corrections to the 
compression modulus B using standard perturbation the- 
ory. The one-loop graphical corrections to B are given 



6B 



d''p[TG\p) + 2A,plG\p)\ pt 



Cd-lQd-l 2 a ( 

{2d -7) ' 



B 



- 7~2d 



(IV.l) 



where is the infrared cutoff of g,,; L will typically 
be of order the system's length in the soft direction, for 
systems that are much shorter in the soft than the hard or 
z directions. For d < 7/2, SB grows with L and becomes 
comparable to S at L = the "nonlinear crossover 
length". This length can be calculated by equating 
\6B\ — B, which gives 



{7 -2d)Bi-f^K-^ 



(IV.2) 



where 6 = r((i/2)r(3 - d/2)/2. The corresponding 
lengths along the h and z axes can be obtained for look- 
ing for the wavenumbers = and = 1/^^ 
at which the B and 7 terms in the propagator G{p) be- 
come comparable to the K term evaluated at This 
gives: = (e)VA. and ^ = {i^f IK: where 
A2 = yJKlB and \h = \jKj^. The conventional har- 
monic theory only applies when at least one of the length 
scales being probed in the soft, hard, or z directions 
is smaller than the corresponding non-linear crossover 
lengths ^ . To go beyond these crossover lengths, we 
employ the RG which has been derived in section Hill 

To do this we use the trajectory integral matching 
formalism [25} ■ as used in(l2j and many other papers. 
In this approach, the RG is used established a connec- 
tion between a correlation function at small wave vec- 
tor in the unrenormalized system, and the same corre- 
lation function at a larger wavevector in the renormal- 
ized system. While the former cannot be evaluated us- 
ing the harmonic approximation theory, the latter can, 
if we rescale to the point where the rescaled wavevector 
lies within the harmonic regime just described. Here, we 
chose the con nected di so rder averaged correlation func- 
tion is C{q) = {\u{q)\^) - {ui-q)) {u{q)) . 

For this correlation function, the connection just de- 
scribed is given by: 



Ciq,B,K,j,At,g,w) 

^[2x+(d-2)ujh+'^z + l]e(j ^ 



qsC ,q^_e ' ,qhe 



B{£),Kii),jii),At{i),gi£),w{i)), 



(IV.3) 



where on the right-hand side the prefactor comes from 
the dimensional and field rescalings. 

First we consider special g 's with q^ = 0, qh = 0. In 
this case we will choose the rescaling parameter is taken 
to he i = £* = In (A/qs), so that the rescaled q^s lies on 
the boundary of the Brillouin zone. Thus, we can accu- 
rately evaluate the correlation function on the right-hand 
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side of Eq. (jIV.3[) using the harmonic approximation, and 
obtain 



C{q,B,K,-f,At,g,w) 

I2x+{d-2)uih+'^z+l]e' _ 



1 



1 



where 



with 



(IV.4) 



(IV.5) 



(IV.6) 



8\/2 27 

In the second equahty of Eq. pV.4p we have ehminated 
£* in favor of qs and used 

Kit) = K^^N^f-^)--+-'+^x-3 ^ 



.{2-d)uj,,-uj,-2x+3--^ 



✓2 



g* / A \ (rf-2)(^h+u;j+2x-3 



(IV.7) 



which is obtained by integrating the flow Eq. (|III.12p . 

The above calculations can be generalized for arbitrary 
direction of q by using a more sophisticated choice of £* , 
in which we rescale until the propagator G{p} takes on the 
same value as its smallest on the Brillouin zone boundary. 
This leads to the condition: 

l{t){qhe^-'')\ (IV.8) 

This choice again ensures that the correlation function 
on the right-hand side of Eq. (|IV.3p can be accurately 
evaluated using the harmonic approximation. Clearly, 
the solution of this equation reduces to £* ~ In (A/qs) for 
the special case qn = O^Qz = 0- In general, the solution 
is 



where we have defined the scaling variables 



and 



X 



Y 



(IV.9) 



(IV.IO) 



(IV.ll) 



The universal anisotropy exponents (z and (^h are given 
by 



a = 2 



(IV.12) 
(IV.13) 



In equation (jIV.9[) . f{X,Y) is a scaling function. Its 
asymptotic behaviors can be deduced by noting that, for 
X ^ 1, F <C 1, it should recover the result for qz ~ 
qh = 0, while for X > l,X^i- > Y^^ t should only 
depend on qz, and for y > 1, Y'^' » X'^'^ , £* should only 
depend on q^- These conditions force f{X,Y) to have 
the asymptotic forms: 

1, X<1, y<i, 

fix,Y) = {x-^, x>i, X?" > r'^^(Iv.l4) 

Y^^, r > 1, Y'^' > 



Using £* as given by Eq. pV.9p and repeating our ear- 
lier calculation for the qh = qz = case, we obtain the 
anomalous elasticity for K and A( for an arbitrary 
direction of q: 



B{q) = B{^^q,y'^{fiX,Y))-^- 

- ieisrfB{x,Y) 



K{q) 



{f{X,Y)f 



{iUsV fK{X,Y) 



Mq) = A, (4%) "{f{X,Y)r 
^ {^qs)-'' ft{X,Y) 



(IV.15) 



(IV.16) 



(IV. 17) 



where we have defined the B, K, and At scaling func- 
tions: 



/B(X,F) = B(/(X,y))- 



fKAX,Y)^{K, At){f (X,Y)r-^^ 
Here, the exponents r]B and r]t arc given by: 



VB 

Vt = 



9i 



9 



32V2 



(IV.18) 
(IV.19) 

(IV.20) 
(IV.21) 



In deriving Eq. pV.17p . we have used the correlation 
function {u{—q)){u{q)) = AtgfG^(g) instead of C{q) , 
since the latter has no dependence on A^. 

In the physical dimension d = 3, to leading order in 
e, r]B = 0.2222, r]K = 0.2963, and rjt = 0.07407 . These 
results are quantitatively reliable, since e = ^, which is 
small. Since both rjB and rjk are positive, the compression 
modulus B and the bend modulus K vanish and diverge, 
respectively, as g — > 0. Physically, this can be understood 
as follows: 

In the presence of layer fluctuations, a compression 
can be relieved by smoothing out these fluctuations, thus 
leading to a vanishing effective compression modulus. 
Similarly, bending a fluctuating smectic leads to com- 
pression of the layers, thus leading to a diverging effective 
bend modulus. 
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V. CRITICAL BEHAVIOR 

In this section we calculate the critical exponents for 
the transition following standard RG procedures [2^ . 
Linearizing the flow Eqs. (|III.19[ IIII.20[ IIIL24I) around 



the fixed point gl 



128\/2 , 
27 



0(e2), D* ^ + 0{e^), we obtain 



dSgi 



M 
dSD 



^^51 
jSg,^ 
Xd6D 



~j6g, 
4:e6g2, 



8V2 



^.91 



8V2 



Sg2 



where 



28 

XD=2--e + 0{e^). 
These recursion relations have the solution: 



^92 

SD 



^^^^ 



36 



16V2 



3KC2 
64^2 



512\/2 

L " ' 

(V.l) 
(V.2) 
(V.3) 

(V.4) 

(V.5) 
(V.6) 

e-^^^ (V.7) 



where the C's are determined by the bare values of the 
coefficients in the Hamiltonian Eq. pil.ip . 

Note that 6D flows to oo for positive C3 > 0, to —00 
for C3 < 0, and to zero for C3 = 0. Thus positive C3 
corresponds to the A phase, negative C3 to the C phase, 
and C3 = to the AC transition. In terms of the bare 
values of 51,2 and Z), the AC transition occurs when C3 — 
0; a little algebra shows that this happens when 



DiTAc) 



3K 
I6V2 



{92 -91): 



(V.8) 



which describes a plane in the gi-g2-D parameter space. 
Unlike in mean field theory, the critical value of Z?, and 
hence the critical temperature Tag depend on both gi 
and 52- 

C3 measures how far the system is away from the 
critical point. Assuming that C3 is a smooth, analytic 
function of temperature near Tag (and we have no rea- 
son to assume otherwise), then near the critical point, 
C3 cx t = This is the sole input we will need in 

our analysis to get the critical dependence of all quanti- 
ties on t. 

For ^ — > 0+, the correlation lengths of the original sys- 
tem and those of the rescaled system are connected by 



6 

6 



(V.9) 
(V.IO) 
(V.ll) 



This allows us to calculate the singular i-dependcnce of 
the correlation lengths. Under the RG, according to flow 



Eqs. dVSlElEZl 5gi 

SD{e) (X te^°^ . At 



0, Sg2 0, and 5D grows as 



A, 



■Int 



(V.12) 



6D{£) is of order one, and the rescaled system becomes 
non-critical with finite (,s,h,z{i)- Plugging £ = £* into the 
above equations we obtain 



with 



1 



C,zVs = 1 - 



7 

27 
7 

27 
10 

27 



e + 0(e2), 
e + 0{e^). 



(V.13) 
(V.14) 
(V.15) 



(V.16) 
(V.17) 
(V.18) 



The temperature dependence of the order parameter 
0jVj the angle between N and z, can be calculated in a 
similar way. For t — > 0^, fl^v is small and approximately 
On = \dsu\. \dsu\\B. the original system and that in the 
rescaled system are connected by 



,(-i+x)f* 



(V.19) 



\dsu{l*)\ can be calculated using mean field theory since 
the rescaled system is deep into the C phase, and given 
by 



\dsu{t)\ 



-2D( 



(V.20) 



Since the dimensionless coupling 2 are finite constants 
near the RG fixed point, and we have chosen our RG 
rescaling factors so that, K{i), B{£), and 'y{£) are kept 
fixed at their bare values, we have, from the definitions 
Eqs. mm mm . w{£*) ^ and g^ ^ 

Furthermore, D{£*) is of order one. Therefore, Eq. 
(rvT20)) leads io \dsu{ £*)\ - ^At(£*). Substituting this 
resul t into E q. ([¥19)) . using the recursion relation equa- 
tion pil.isp for At; and using equation (|V.12I) for £*, we 
get 

\dsU\ |t|^[(2-f')Ch-C=-7+4d-?)i5-T,t-(d-3)?)A-] (V.21) 

which implies 

p = 4[4d-10 + (5-d)?7K-r;B-2ryt] 



1 
1 2 

2^9 



e + 0{e^). 



(V.22) 



The first and second equalities are respectively the exact 
scaling relation p.l7p and the e-expansion result (|I.2[) 
quoted in the introduction. 
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The specific heat exponent can be obtained by similar 
techniques |25|. The disorder averaged free energy of the 
system is invariant under the RG, that is, 

F = Vr7r[K{n,B{n,iin,D{t),Mn] 

+AF{e*), (V.23) 

where Vr and are the volume and disorder averaged 
free energy density of the rescaled system, respectively, 
and AF{£*) represents the contribution to the free energy 
from the degrees of freedom already integrated out in the 
RG process. We will focus in what follows on the first 
term on the right hand side of equation (|V.23p . since it 
scales the same way as the full free energy. If renormalize 
for the same renormalization group time £* that we used 
above to calculate the correlation lengths, we will again 
have renormalized out of the critical region, which will 
permit us to calculate the free energy in mean field the- 
ory. With this choice, Vr is related to the volume of the 
original system by Vr = T/oe"[(''"2)"''+"-+il''* . Further- 
more, because we have renormalized out of the critical 
regime, fr only depends on K{e*), B{t), 7(r), £»(£*), 
and At(£*) since it can be evaluated by ignoring the an- 
harmonic terms. The only singular temperature depen- 
dence oiJrin comes from At(r), since K{t), B{£*), 
and 7(i!*) are fixed at their bare values and D{£*) is of 
order one. This singular temperature dependence can be 
calculated by noting that is linear in At{£*); that is 



dAt{£*) 



n->-0 (2'!T)"-n J ^ — ' 
a/3 

(V.24) 



constant. 



Therefore, 

F cx Vbe-[('^-2)'^''+"=+il'^' [At(r) +cons.] 

where in the second line we have used the solution of 
the recursion relation (jIII.15|) and the expression (jV.12[) 
relating £* to the reduced temperature t. 

Because the specific heat is the second temperature 
derivative of the free energy, the power of t on the right 
hand side of the bottom line of equation (|V.25P is just 
2 — a, where a is the usual specific heat critical exponent. 
Thus, 



a = 2 - ^ [4d - 6 + (3 - d)m - 2r,t - vb] 



(V.26) 



The first and second equalities are respectively the exact 
scaling relation (jl.lSp and the e-expansion result p.l6[) 
quoted in the introduction. 

Once we probe the system at length scales longer than 
the correlation lengths ^s,h,z in the corresponding direc- 
tions, the elastic constants and disorder variances become 



wavevector independent. This can be seen either by re- 
peating the RG analysis just done above for small D for 
large D; alternatively, one could repeat the perturba- 
tion calculation equation (|IV.ip with a nonzero D, and 
find that all perturbativc corrections now converge in the 
infra-red. 

Instead of wavevector dependent, however, all of the 
elastic constants and disorder variances now become the 
temperature dependent. We can obtain this temperature 
dependence by matching. 

For example, the temperature dependence of B can be 
obtained by matching, right ai Qs = , the wavevector 
dependence oi B{qs,qz = 0, (f/i = 0) for qs » £,J^ given 
in equation (jIV.15|) onto the wavevector independent be- 
havior we expect for qs <C ^7^- Doing so, we obtain 

BlT) = B{if/UTW^- 

Wc can obtain D{T) using the definition of the corre- 
lation length ^s, which implies K{q)q'^ = D{T) at qs = 
C\q. ^0,qh = 0; this implies DiT) (x (Cf 

Similar arguments give K{T), and the disorder vari- 
ances At_c(T), where the needed critical renormalization 
of the compression disorder variance Ac{q) is computed 
in appendix (|B|). Summarizing all of these results: 



B{T) = «(r-r^cr 



-11 K 



« {T-TacY 



(V.27) 



(V.28) 



D{T) = D{i^ lisf-'^"-^ ^{T-TacT-^'^'^''^ , 

(V.29) 



AtiT) = At(Cf/6)"'"«(r-T^c)-''*'^' 



A,(T) = A,(Cf/$.)-''=cx(r-T^c)-''=''= 



(V.30) 



(V.31) 



These temperature dependences are needed to predict 
the singular dependence of the light scattering in the A 
and C phases near the critical point, as will be shown in 
section (|VII|) . 



VI. PROPERTIES OF THE A AND C PHASES 

The analysis of the preceding two sections has focused 
on the critical region near the AC transition. What about 
the properties of the A and C phases themselves in this 
random biaxial environment? 

These prove to be quite interesting and novel as well: 
we find that both the A and C phases behave, at long 
distances, like an XY ferromagnet in a random field jT8|. 
This should be contrasted to the behavior of this system 
in uniaxial disordered environments, in which case the A 
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phase also looks like an XY model in a random field, but 
the C phase belongs to the far more disordered "m=l 
Bragg Glass" universality class [l7j. 

We begin with the A phase. At long length scales 
rs < Cs, rz < £,z, and rh <. ^h)(i-e., D > 0), the effective 
Hamiltonian for this phase is: 



H = d"^ 



+Vp{u^cl,{f))]. 



(vi.i) 



Note that we have dropped the bend term be- 
cause it is less important than the (dsu)^; near the tran- 
sition, we could not do this, since D{T) vanished there. 
But in the A phase, since D > 0, we have this term, 
and it dominates the bend term at long length scales as 
defined above. 

An additional effect of this extra term is that it sup- 
presses fluctuations of u sufficiently that the g and w 
anharmonicities in (jll.l[) are also irrelevant, as are the 
random tilt terms. Hence, the only relevant anharmonic- 
ity is that implicitly contained in the random pinning 
potential Vp{u ~ '^(f')), which we remind the reader is a 
periodic function with period a. 

Note that neither the compression modulus B{T) nor 
the coefficient D{T) in the Hamiltonian Eq. (jVI.l[) are 
given by their bare values. Instead, both are renormal- 
ized by the critical fluctuations and become non-trivially 
temperature dependent. They are given by Eqs. (jV.27[ 

EMI). 

Like the model for the smectic A phase in uniaxiallvfv^ 
stretched aerogel, this model is effectively a random field 
XY model (RFXY), with an anisotropic spin wave stiff- 
ness. 

To show this, we rescale lengths as follows: 



and replace u by a rescaled variable 

6{v') = q^u(r) , 



(VI.2) 



(VI.3) 



where go = 27r/a is the reciprocal lattice constant. This 
latter transformation has the property that the physical 
state of the system is now invariant under the transfor- 
mation: 



e ^ e + 27rn 



(VI.4) 



where n is any integer, since the smectic state is invari- 
ant under translations by any integral number n of layer 
spacings a. This is precisely the symmetry obeyed by 
the angle variable dxY giving the angle between the lo- 
cal spins and some reference direction in a ferromagnetic 
XY model. And, indeed, in the rescaled variables, the 
model becomes: 



K 



XY 



V'6 



XY I 



- (t^XYif)) 



(VI.5) 



^ffhereV^^ {6) = Vp{^) is a periodic function with pe- 
riod 27r, (pxYi'r) = ^^^'•''•^ is a short-range correlated ran- 
dom variable uniformly distributed between and 27r, 
and, in spatial dimension d = 3, the spin wave stiffness 

KxY = 

This Hamiltonian (jVI.SP is easily recognized |26| as just 
the long- wavelength Hamiltonian for a ferromagnetic XY 
model in a random field, which has been extensively 
studied (3]. Fisher [l^ showed that , in this model, the 
spatio-temporally Fourier transformed two point corre- 
lations of are given by {9{q')6{—q')) = ^pz^, where 
C{d) is an constant of order one and C(3) ~ I.IOtt^. 
Using the rescaling (jVL3[) with this result then implies: 
{u{q')u{-q')) = Then using u{q') = ^W/Vu{q), 

where V and V are respectively the volumes of the sys- 
tem in the original and primed coordinate systems, and 
rescaling q' back to q using the inverse of the rescalings 
(|VL2l) . we find, in d = 3, 



{u{q)u{-q)) 



[Bql + Dql+WhY ql 
The real space correlation function 



C{r) = {[u{r) - u(0)]2) 



(VI.6) 



(VI.7) 




can be calculated using a similar strategy. The real space 
correlation function governed by Hamiltonian (jVI.5|) has 
also been derived by Fisher |18|. and is, in d = 3, 

{[e{r') -6l(d)]2) = 1.101n(r76'), where h' is the cutoff 
length in r ' space. Then rescaling 9 and r ' back to u and 
f , respectively, using the rescalings described above, we 
obtain 



(VI.. 



where &<,.h,z are respectively cutoffs along s. h, and z in 
r space, and related to b' by the rescalings (jVI.2[) . which 
give b, = b'^D/B, bh = b'yfijB, and bz = h' . 

The physically more important correlation is Fn{r) = 
^gjnqo[u(r)-u(o)]^^ which can be experimentally measured 
by X-ray scattering, as will be discussed in next section. 
In the approximation that the fluctuations of u are Gaus- 
sian, which Fisher [isj has shown to give the correct power 
law decay of the XY analog of Fn{r) to leading order in 
e = 4 - d, 



Fn{r) « cxp 



9 9 



C{r) 



(VI.9) 



Plugging Eq. fVLS)) into Eq. (|VL9| we get 



2 / \ 2 / N 2 ' 




-0.55n^ 



(VI.IO) 



14 



i.e., quasi- long-ranged translational order. However, un- 
like in bulk sniecticsj2^. the power law exponent has no 
dependence on either the temperature or the quantitative 
elastic properties (e.g., the bare values of the elastic coef- 
ficients), and is therefore universal. Furthermore, Fn{r) 
is isotropic in scaling, unlike the clean casej2^. 

Now we discuss {u{q)u{—q)) at small g's in the C 
phase. Since £> < 0, u = is not the disorder aver- 
aged ground state. Therefore, we need to first expand 
the Hamiltonian around the state dzU — — „ ;> , dhu = 

0, d.u -- 



DB 



Bw-g' 

which is one of the disordered aver- 



aged ground states in the C phase, as obtained by min- 
imizing the Hamiltonian Eq. pi.ll[) with the disordered 
terms excluded. We will denote the displacement of the 
layer from the ground state position as u'(f). Expanding 
in small u ', and ignoring all irrelevant terms, we obtain a 
Hamiltonian which is essentially the same as Eq. (jVI.l[) 
but with u{r) replaced by w'(r). Thus, we conclude that 
the universality class of the C phase is als o that of the 
random field XY model, and {u'{q)u'{—q)) is also given 
by equation (|VI.6p . Likewise, real space correlations are 
also given by equation (|VI.10[) . 



VII. PREDICTIONS FOR SCATTERING 
EXPERIMENTS 

In this section we work out the implications of our 
theory for both visible light and X-ray scattering exper- 
iments. 



A. X-ray Scattering 

X-ray scattering experiments provide a direct experi- 
mental measure of the density-density correlation func- 
tion {p[r)p(f')) , where p{r) is the molecular number den- 
sity. This correlation function is an important quantity 
in smectics, since it reflects the translational order of the 
system. The X-ray scattering intensity is related to its 
Fourier transform by 



I{q) cx / d'^r'd'^r" {p{r')p{r"))e-'^<^'' 



(VH.l) 



In a smectic, p{r) can be expanded in a Fourier series 
with period a, the layer spacing between nearest layers. 



(VH.2) 



where p„ is the complex amplitude of the nth harmonic 
of the density wave. In writing this equation (|VII.2p . 
we have included the effect of fluctuations u(r) in the 
positions of the layers, but not those of the magnitudes. 
The latter fluctuations can be shown [27| to be irrelevant, 
in the RG sense, close to the transition. 



Inserting Eq. (|VII.2p into Eq. (|VII.1I) we obtain 

I{q) oc ^ |p„|2 f d\ e*("«°^-«^-^X(f), (VH.3) 

where 

Fn{r) = (e™'3o[«(*=')-«(o)]). (VII.4) 



In deriving Eq. (| VII. 31) . we have changed variables of 
integration from r ' to r ^ r' — r" and used the fact 
that, for homogeneous systems. 



l^^inqo[u(r')-u{r")]\j _ ^gingo ["(r) -u(0)] ^ ^ (VII. 5) 

Because all of the Fn{r) are slowly varying functions of 
position (because u{f) itself is), the sum in Eq. (|VII.3[) 
is dominated by the term with nq^z the nearest Bragg 
peak to q. Therefore, the X-ray scattering intensity near 
the nth Bragg peak can be simply written as 

In{5q) oc \pnf I d^r e'*«-^X(f), (VH.G) 



where 5q is defined as the deviation from the center of 
the nth Bragg peak, 5q = g — nq^z. 

We have already obtained the scaling of the limiting, 
large r behavior of F„ (r) in both the A and the C phases 
in section (|VI|) . 

To get an idea of the behavior of Fn{r) for smaller r^s, 
which proves to be extremely important for getting the 
full scattering, we will begin by approximating Fn{r) by 



cxp 



? 9 



C{r] 



(VH.7) 



where 



C{r) {[u{r)-u{0)f) 
2 



(27r) 



^ , dU [1 - cos [q- r)] {u{-q)u{q)). 



(VII.8 



This expression would be exact if the fluctuations of the 
displacement field u were Gaussian. As discussed above, 
this approximation can be shown (isj to give the correct 
scaling of F„ (r) at the longest length scales in both the 
A and the C phase. Its validity for length scales shorter 
than the correlation lengths (i.e., for «C ^s, rz <C ^z, 
and r/i <C S,h) is more questionable; here we will begin 
by assuming equation (jVH.Sp and exploring the conse- 
quences of that assumption. Then we will discuss possi- 
ble perils of the approximation (jVII.8[) . and argue that 
the qualitative form of the scattering is given correctly 
by it, although certain quantitative results (in particular, 
the precise numerical value of the universal exponent 
in equation p.23|) for the spike width in the X-ray scat- 
tering), may be incorrectly given by it. 
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We can calculate C(r) in the Gaussian approximation 
using our knowledge of the disordered averaged fluctua- 
tions in momentum space, {u{—q)u{q)). 

At large qs (i.e., qs > C\ or qz > or > 

^^^), the term (dgu)'^ in the Hamiltonian Eq. (jll.lip is 
subdominant to and hence negligible. In this case 

the random pinning disorder is irrelevant, and we treat 
the tilt-only model using the RG derived in Sec. IIIII 
{u{—q)u{q)) can be calculated by following the method 
described in Sec. ITVland given by 



{u{-q)u{q)) 



[B{q)q'z 



K{q)qi 



(VII.9) 



where the wavevector dependences of B{cj), K{ct), and 
At{q) are given by equations (|IV.15IIV.16IIV.17| . respec- 
tively. 

Plugging {u[—q)u{q)) into Eq. (|VII.8|) we calculate 
C(r) asymptotically. For r less than the correlation 
lengths (i.e., rs.h^z ^ S,s,h,z), the dominant contribution 
to C(r) in Eq. (|VII.8p comes from the integral over large 
(f's (i.e., qs > C\ or qz > C\ or qn > C^), and C(r) 
is thus given by 



C{r) « Xl 



[(If) 



r/Ch 



i^)^ '(if) 
W ^ iw) ' (#) 



(VII.IO) 



where F = 1 -I- 774 — + ^ . Plugging this result into Eq. 
(jVII.7p gives an exponentially decaying translational cor- 
relation function; i.e., short-ranged smcctic translational 
order. 

For r larger than the correlation lengths (i.e., rg 3> ■fs, 
or > or rh > 61), C{r) crossovers to 



C(r) 




(VII.ll) 



The first piece comes from the integral over small g 's 
(where we have used "small" in the sense described ear- 
lier). It follows from Eq. (|VI.8p with the cutoff lengths 
bs,h,z replaced by correlation lengths S,s,h.z, respectively. 
The second piece comes from the integral over large 
g's and connects with C(r) Ea. (|VII.10[) at the crossover 
point. 

The x-ray scattering intensity near the nth Bragg peak 
In{Sq) can be calculated using Eq. (jVII.6|) . Plugging the 
above calculated C{r) into Eq. (jVII.7|) . we find F„(r) 
has the following behavior. For r less than the corre- 
lation lengths Fn{r) decays exponentially fast, while for 
r larger than the correlation lengths Fn (r) decays much 
slower and as a power law. Thus, in the large Sq limit 



the main contribution to In{5q) Eq. (|VII.6p comes from 
the integral over small r. This implies that far away from 
the center of the Bragg peak, the line shape of the scat- 
tering intensity is qualitatively Lorentzian-squared. The 
line widths of this squared-Lorentzian are anisotropic and 
given respectively by 



= (c)-'-(e)"'(^) ^ 



-2G,„/r 



(VII.12) 



(VII.13) 



where ^ are the line widths of the correlation function 

Fn{r) and obtained via C(r) = a? . In the limit 5q — >■ 0, 
IniSq) converges for n > 2 and diverges for n < 2. There- 
fore, only the first two Bragg peaks diverge, following the 
power law: 



IniSq)(XZUn (Sqs^s) +iSqh^h) + (SqzCzY 



where n7„ = S,sS,z^h exp 



(VII. 14) 



has a ex- 



in'>^lqUUe) 
ponentially strong dependence on the temperature, which 
is responsible for the fast vanishing of the divergent 
Bragg peak approaching the critical point. The diver- 
gence shown by (|VII.14p comes from the integral in Eq. 
(|VII.6P over large r . 

The crossover between these two distinct scattering 
patterns can be estimated as the follows. Since the di- 
vergent (power law) Bragg peak vanishes extremely fast 
approaching the critical point, we expect that near the 
critical point, the crossover occurs at very small Sqs which 
are much less than the widths of the Lorentzian square. 
Thus, at this crossover the Lorentzian part of /„((Sgc) is 
approximately proportional to £, f£,:^f- Equating this to 
the power law part given by Eq. (jVII.14p we get, for the 
nth Bragg peak. 



Ss 'ih Sz 



137,, 



exp I 



Arit 



(VII.15) 



where An is non-universal, VIg (where the subscript 
"G" stands for "Gaussian" ) and T" ^ ^ are universal and 
given, respectively, by 



Tn 



1 
2 



+ —e + 0{e^), (VII.16) 



27 



Vs^hyz 
-3+ .55n2' 



(VII.17) 



These results were derived in the Gaussian approxi- 
mation. Since we have shown that the smectic elasticity 
and fluctuations are controlled by a non- Gaussian fixed 
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point, one might well question the vahdity of the Gaus- 
sian approximation. 

How could we go beyond the Gaussian approximation? 
In principle, this eould be done by integrating all of the 
degrees of freedom out of the Hamiltonian (completeH) 
except for the single variable that enters the X-ray scat- 
tering near the n'th Bragg peak, namely: 



Xn{r) = nqo[u{f) - u(0)] . 



(VII.18) 



The hypothetical integration out of all other variables 
would leave us with a known probability distribution 
P{xn{r)) of this single variable a;„(r). 

In reality, of course, it is far too difficult to perform this 
hypothetical integration. We do, however, know some- 
thing about P{xn{r)): its variance, which is just C(r), 
equation (jVII.8[) . This sets the scale of a;„(r) fluctua- 
tions; hence, we can write P{xn{r)) in a scaling form: 



P{x^{r)) = Vke 



where 



fc(r) 



0(1) 



{nqorCir) 



(VII. 19) 



(VII.20) 



If we take U{y) = y^, then we recover the Gaussian 
theory just presented. More generally, we can write: 



/oo 
dx 
-OO 



(VII.21) 



For |r| large, A;(r) is small, and one can accurately eval- 
uate this integral, for a given U (y) , using the method of 
saddle points. That is, one finds the place where the 
argument of the exponential in (|VII.21[) is extremized, 
expands about that point, and performs a Gaussian in- 
tegral using this expansion. In general, locating the sad- 
dle point requires analytically continuing the integral in 
(|VII.2ip into the complex plane. 

The first step is, of course, locating the saddle point. 
Defining 



f{x) = ix~ , (VII.22) 

this point is where 

,f'{x) = i - ^/k{¥)U'{^/k(f}x) = . (VII.23) 

For example, consider U{y) = y'^jl^. The condition 
(|VII.23|) becomes 



whose solution is [2£ 



x± 



(32fc2)3 



(VII.24) 



(VII.25) 



To exploit these saddle points, we deform our contour of 
integration from the real axis to a line parallel to the real 



(-3'°1) 


y 

(3'°,1) 
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FIG. 6: The complex contour we choose for the saddle point 
evaluation of the integral ()VII.2H) for the translational corre- 
lation function Fn(f). 



axis, but displaced up (i.e., in the imaginary direction) 
from it by the imaginary part of x± , as illustrated in Fig. 
[51 (The contour is connected to the original contour by 
vertical segments at ±c» that contribute nothing to the 
integral.) This passes the contour of integration through 
both saddle points x±. 

Since the integral is now dominated by the region near 
these two saddle points (as we will verify a posteriori in 
a moment), we can expand the integrand around each of 
the saddle points. Noting that 



and 



-Uk^xl 



(1± V3i)(VII.26) 



ix± 



i{±V3i ~ 1) 
8(2fc)i 



we can use the expansion 



f{^) « fix±) + -f"{x±W 



(VII.27) 



(VII.28) 



with u real, near each of the saddle points. Making this 
approximation, and extending the region of integration 
over u near each saddle point to run from — oo to oo 
(an approximation that we will justify for small A: in a 
moment), we are then left with two Gaussian integrals, 
whose sum is Fn{r). This gives 



Fnir) 



exp 



8(2fc)i 



3^3 
8(2fc)i 



(VII.29) 



Now, to justify the saddle point approximation: 
This approximation will be valid provided that the ex- 
pansion (|VII.28|) is valid over the range of u that dom- 
inates the integral near each of the saddle points. This 
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implies that it must work for 

The neglected next order term in the expansion (|VII.28P 
is ~ f'"{x±)u^ ^ — — ^-^'^l , where we have used the esti- 

mate (jVII.30p for the range of u that dominates the in- 
tegral near each saddle point. As long as this is <C 1, 
the saddle point approximation will be good. Since, 
from (jVII.26|) . f"ix±) - fci while f"{x±) ~ k^x± ~ 
k'^k~3 fc3 (where wc have used (|VII.25P for x±), we 
f"'ix±) „, ki 

k 



with 



have 



ks <C 1, provided fc <C 1, as it is 



/"(^±)^ 

for large as noted earlier. 

Thus, for large |r|, the saddle point approximation is 
valid. 

More generally, if we take U{y) = y", then we get, by 
exactly the same sort of steepest descent analysis, 

i^„(r)occxp(-0(l)xfc-^(")) (VII.31) 



with 



e(n) 



2(n-l) 



(VII.32) 



This is readily seen to reproduce the results for the 
Gaussian case (n = 2), and the quartic case we just ex- 
amined in detail (n = 4). Note that for all values of 
n > 2, the exponent O is bounded above and below: 



(vn.33) 



with the upper bound satisfied by the Gaussian case, 
while the lower bound is approached in the limit n — >■ oo. 
Since we expect n > 2, since H contains higher powers of 
u and the value of Q is determined by the highest power, 
we expect Q to satisfy (|VII.33p . 

As a final example, let us consider U{y) — e^^. The 
steepest descent analysis now gives 

Fn{r) cx exp (^-0{lnk) x fc"^) , (Vn.34) 

which is effectively (up to slowly varying logarithmic cor- 
rections) of the form ()VII.3H) with Q — ^] i.e., just the 
limiting value of Q{n) as n — >■ oo (as is hardly surprising, 
since the Taylor series for involves arbitrarily large 
powers of y). 

While hardly exhaustive, these examples give us confi- 
dence that the leading behavior of F„ (r) is always given 
by (|VII.3ip . with 6 satisfying (|VII.33|) . Using this, and 
the definition (jVII.20|) . and repeating the analysis given 
above of the Gaussian case, we find that (S<?"^ ^ is given 
by: 



CX t' 



Ss ?h Sz 



cxp {-A„t^^) 



(VII.36) 



where flc is give n by (|Vn.l6|) . T by (|Vn.l7|) . and 9 
obeys the bound (jVn.33|) . Combining this bound on Q 
with an estimate of the uncertainties in the e-expansion 
result (jVII.16p for flo leads to the bounds on ^l quoted 
in the introduction. 

Near the critical point the dominant i-dependence of 
Sq^l 2 is exponential. Thus, a plot of lnln(5g"'^ ^ versus 
Ini should give a straight line with slope fl. 



B. Visible Light Scattering 

The fiuctuations of the nematic director scatter the 
light strongly. The scattering intensity is proportional to 
the fluctuations [2411: 



I{q) (X (ei-fto) (\Sn{q) ■ erl'^) + {er ■ no) {\Sn{ci) ■ ee\'^) 
+2 (e£ • no) {er ■ no) ■ {Sn{q)Sn{-q)) ■ er, 

(VII.37) 

where ei and Ct are respectively the polarization direc- 
tion of the incident and the transmitted light, and Sn is 
the deviation of the nematic director from its averaged 
value no: Sn = n — ho- 

In the A phase no = ^, so the largest components of 6n 
are those d7i± perpendicular to z. Thus, the scattering 
intensity is proportional to a linear sum of {5nj^5n^)] 
namely. 



m ^ J2 D..j{5ni{q)5nj{-q)), (¥1138) 

where the coefficients Dij are of order one and solely de- 
termined by the polarization directions of the incident 
and transmitted light, and 5nj^ denotes the i 'th compo- 
nent of 5n^. 

Since, in Sec. [ill we have shown that the fluctuations of 
the linear combination Sn' ± (equation pi.Sp ) of 5nj- and 
(9j^u are massive, {Snl{-q)5nj-{q)) w {d^u{-q)d--u{q)) , 
where d^- denotes the partial derivative along the i 'th 
component in the _L plane. Therefore, equation (|Vn.38p 
can be rewritten as 



where 



a,iq}^L,,{q)qi{\uiq)\^^ 



(vn.39) 



(VII.40) 



(vn.35) 



with Lij{q) = q^qj-fq"^ the projection operator along _L, 
and the coefficients Eij are independent of q and of order 
one. 

An expression for (|u((7)P) has been obtained in Sec. 
IVII Al However, those results have to be improved to pre- 
dict the light scattering. Specifically, in addition to the 
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fluctuations induced by the random tilt in the s direc- 
tion and the random field, we need to include in 
the fluctuations induced by the random tilt in h direc- 
tion and the random compression, as embodied in the 
and Ac terms of equation (jIL12[) . These last two types 
of disorder, though " irrelevant" in the RG sense, have 
experimentally detectable effects on light scattering. 

If any component of q is large compared to the corre- 
sponding inverse correlation length, (i.e., if ^ ^7^, or 
Qh 3> or Qz 3> ^J^), but all components are small 
compared t o the ultraviolet cutoff A, we find that the 
< |u((7)p > in (|VII.40|) is given to a good approximation 
by 



< Hq}\^ > = 



At jq) ql + Ahql + [q) ql 
{B{q)ql+^ql + K{q)q^f 



(VII.41) 



where the g-dependences of B{q), K{q), and At{q) were 
derived in section IIV| and are given by Eqs. ()IV.15| 
IIV.16[ and IIV.17[ respectively. The wavevector depen- 
dence Ac{cj) of the z component of the random field h 
is derived in Appendix B, and given by equation (|B.7I) . 
The correlations Afi{q) of the h component of the ran- 
dom field h do not acquire any anomalous wavevector de- 
pendence since there is no graphical correction to Ah (q) ■ 
The fluctuations caused by the random field disorder are 
not included in Eq. (jVII.4ip , since their contributions to 
< |w((7)p > arc always dominated by the others for large 

q- 

Likewise, if all of the components of q are small com- 
pared to the corresponding inverse correlation length, 
(i.e., if qs < and qu < and < ^J^), we 

obtain 



< Hq)? > 



At {T)ql + Ahql + A, {T)ql 



1.10TT^y^B{T)D{T)j 



{B{T)ql + -fql + D[T)qlY {B{T)ql + 7?^ + D{T)ql)-- ql 



(VII.42) 



r 



where the third piece comes from the random field disor- 
der. Recafl thatB(r), K{T), D{T), At,c{T) arc no longer 
wavevector dependent in this regime of wavevector, , but 
are renormalized by the critical fluctuations, and, hence, 
temperature-dependent; this temperature dependence is 
summarized in equations lv:27llv:^ lv:29l Iv30l IV^Tl 

As mentioned in the introduction, These rather compli- 
cated expressions (equations (jVII.4ip and (|VII.42I) ) sim- 
plify considerably in a variety of natural limits. 

Consider first varying wavenumber q for fixed, generic 
direction q of the scattering wavevector q. By "generic 
direction q " , we mean a direction of q for which no com- 
ponent of (7 is ^ 1. 

Consider first the regime of s large compared to the 
correlation lengths, (i.e., qs ^ S,"-^, or q^ » or 
1z ^ but small compared to the ultraviolet cutoff 

A, in which equation (jVII.41[) is valid. Since all three 
components of q arc comparable for a generic direction 
of q, this amounts to saying that the magnitude q of q 
must be large compared to the smallest of the inverse 
correlation lengths qt ^ (.^"^ , and qz ^ Cj^: which 
means it must be large compared to the inverse of the 
biggest correlation length, which is ^/j, since ^/i cx ^l*" 3> 
^2 oc ^5"^ I where the inequality follows from the facts that 
(^ = 2- ^>G = 2- (2ii+!ZB) and ^ 00 as the 
transition is approached. 

Thus, we arc. for a generic direction of q, in the regime 
in which equation (jVII.41[) is valid, provided 



(VII.43) 



where Vh is given by equation p. 6 



For q 's satisfying the condition (jVII.43p , with a generic 
direction of q, the dominant cutoff in all of the scaling 
functions for all of the anomalous elastic coefficients in 
equation (jVII.41[) is qh. 

To see this, begin by noting that the scaling arguments 



X 



N 



and 



Y 



N 



fee) 
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(VII.44) 



(VII.45) 



that appear in all of the scaling functions (e.g., equation 

00 and 



(jrO]) for e*) are both > 1, since X oc q^'^' 
Y oc q^^'''^ — > 00 as g — > 00, since both C,z and <^h are 
> 1. This implies that either qh or qz dominates the 
scaling functions. To determine which, we form the only 
possible gs-independent ratio of X and a power of Y: 



X 

Y 



oc q 



as g 



The 



vanishing of this ratio implies that qh is the dominant 
cutoff in all of the scaling functions. 

Using this fact implies that, for the generic directions 

of q that we are considering, At(g) oc g ^ and Ac(g) oc 
— ^ 

q '^h , This in turn implies that the Ac(g) term in the 
numerator of eqn. (jVII.41| dominates the At{q) and Ah 
terms, since the ratio of the At{q) term to the Ac(g) 

term vanishes like g , which vanishes as g — )■ 0, since 
rjc — rjt — 2 — rjB — rjx > 0; likewise, the ratio of the A^ 
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term to the ^c{<l) term vanishes hke which vanishes 
as <; — >■ 0, since rjc > 0. 

Similar arguments show that the 7 term, which scales 
like for the generic directions of q that we are consid- 
ering, dominates the denominator of eqn. (jVII.4ip . 

Taking these two terms as dominant in the numerator 
and denominator of eqn. (jVII.41[) . respectively, we find; 



< \u{q)\^ > cx g-^-TT-^ , q » , (VII.46) 



.4+ 
1h 



for generic directions of q. Using this result in equation 
(|VII.40p . we obtain for the light scattering intensity, in 
this regime of wavevector: 



Iiq)(xq 



4+^ 



1 » ' 



(VII.47) 



1h 



Using the epsilon-expansion results (|I.13|) and (|I.15[) for 
rjc and (h, we obtain for the e-expansion for the power 
law in this expression: 



Hi 



8e 
27 



(VII.48) 



As discussed above, once qh decreases to the cor- 
relation function < |ti((7)p > crosses over to the form 
equation (|VII.42p . In this regime, as mentioned earlier, 
all of the elastic moduli and disorder variances cease to 
be wavevector dependent, but become temperature de- 
pendent, as given by equations (|V.27IIV.3ip . As q first 
decreases below the Ac term continues to dominate 
< >; in this regime, this implies < |u((7)p >(x ^, 

and, hence, that the light scattering intensity I{q) is con- 
stant. Eventually, however, as q continues to decrease, 
the "random field" term (i.e., the y/TTW) term, must 
start to dominate, since it diverges like ^ as g — 0. 
For smaller g 's, < >oc ^, and, hence, that the 



light scattering intensity I{q) obeys 
I{q) (X . 



(VII.49) 



We can estimate the critical value q^ of q at which this 
crossover between constant I{q) and I{q) cx q~^takes 
place by equating the Ac term and the random field term; 
i.e., by setting: 



Ac (T) (gf )2 _ ^B{T)D{T)^ 



F\2\ 



(VII.50) 



where we have used the fact that the q^ term continues to 
dominate the denominators of all terms for this generic 
direction of q, since its coefficient 7 is the only one in the 
denominator that is not renormalized to small values by 
the critical fluctuations. Solving equation (|VII.50p for qc 



gives: 



^B{T)D{T)^ 
{/^c{T)mi)<i 



cx 



\T-Tac\'^ , (VII.51) 



where the universal exponent </> is given by 



^IB - Vk 
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(VII.52) 



where the first equality is obtained using the known 
temperature dependences equations (|V.27p . (|V.29p . and 
(fOI|) for B{T), D{T), and Ac(T), while the s econd 
follows from the known e-cxpansions p.l0[) . (jl.ll[) . and 
p.l3|) . for riK,r]B, and rjc- Equation (jVII.52[) is, of course, 
precisely equation p.27p of the Introduction. 

To summarize, for generic directions q of q, the light 
scattering is given by equation p.25[) of the Introduction, 
and plotted there in Fig. [3] 

As mentioned in the introduction, more information 
can be obtained by restricting the scattering to the sub- 
space of q with qh ~ 0. In this case the only nonzero 
Cij{q) is Css{q)- Therefore, equation (jVII.39|) implies 
that 



I{q)^Css^qi{\u{q}\^^ 



(VII.53) 



As we just did for the case of a generic direction of 
q, here, too, we can distinguish between the cases of tf' s 
large compared to the correlation lengths, (which now 
means qs ^ ^7^, or ^ ^7^)^ but small compared to 
the ultraviolet cutoff A, and th e reverse . In the large 
q case, the correlation function {\u{q}\'^) takes a simple 
scaling form: 



(VII.54) 



where X is the scaling ratio defined in equation pV.lOp . 
This can be seen by noting that, when q^ ~ 0, both the 
numerator and denominator of (|VII.4ip take a scaling 
form. The numerator is 



At (q) qt + Ac (g) qi = /t(X) + q^'^'qtMX) 



iftiX) 



(VII.55) 



where we have used the scaling forms pV.17|) and (|B.7I 
for At and Ac. 

Using the exact scaling relation p.l9[) for 77c, we find 
that 77t - 2 - ?7c = -4 + ?7B + tjk = -2 (2 - 
—2(^2, where we have also used the exact scaling law C 

9 _ Vb+Vk 
^ 2 

Using this in (|VII.55[) implies that 



Atiq)ql + Ac(g)g; 



Mx)+^ Ux) 



— 'ds J a 



AX) 



(VII.56) 



where we have defined the scaling function for the nu- 
merator 



fnnrniX) = MX) + 



(e 



-I 2 



X^UX) ,(VII.57) 
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and used the definition (jIV.10[) of the scahng variable X. 

Similarly, the square root of the denominator of 
()VIL4ip can be written, when qh = 0, as: 



K{q)q. 



(VII.58) 



where we have used used the scaling forms (|IV.15j) and 
PV.16P for B and if, we have defined the scaling function 
for the denominator 



/de„(X) = /x(X) + 



-I 2 



N 



X\fB{X) (VII.59) 



and we have once again used both the definition (|IV.10p 
of the scaling variable X , and the exact scaling relation 



2 _ Vb+Vk _ 



Using the final equalities of (jVII.55[) and (jVII.58[) in 
(|VII.4ip implies the scaling form equation (|VII.54p , with 



fu{X) 



fn\im{X) 



(VII.60) 



and used the definition (|IV.10p of the scaling variable X. 

The scaling function fu{X) must have limits such that 
the full (|u((7)p) correlation function becomes a function 
only of Qs for X 1, and only of for X ^ 1. This 
implies, by arguments virtually identical to those used 
in, e.g., section (jIV[) . 



I constant, 

A ^2 



x< 1 

X > 1 



(VII.61) 



which in turn implies 



9; 



'(VII.62) 



Using equation (|VII.53p . we can obtain the light scat- 
tering intensity from this: 

l[q) (X \ 2„fc-6-,jt (Vll.oS) 



qhz 



crossovers as we decrease q: first, at g^.z = ?7z, at which 
point all of the parameters c, -B, and K cease to be 
wavevector dependent; (i.e., when the crossover between 
Eqs. (|VII.4ip and (|VII.42p occurs). Then, as we continue 
lowering q, at a smaller a crossover from the u — u cor- 
relations being dominated by random compressions and 
tilts to being dominated by the random field term (i.e., 
crossover from the first to the second term dominating in 
equation (|VII.42p . Taking all this into account leads to 
different behaviors of (|m((7)P) in the six different regions 
of the qii-qz, qh = plane illustrated in Fig. 21 
The locus FG, given by 



FG: qzi 



N 



(qsi 



(VII.64) 



is where the crossover between the two limits of equation 
(|VII.63P occurs. 

The locus EF, given by 



EF: gz-C\ 



(VII.65) 



is where we crossover between expressions Eqs. (jVII.4ip 

and (IVII.42I1 on the g^Cf > (^s^f side of FG occurs, 
while DF, given by 



DF: qs=Cs 



(VII.66) 



plays the same role on the q^S,^ ^ {q.sS,^)'^' side. 
The locus OF, given by 



OF: qz = ^qs, 

S2 



(VII.67) 
(VII.68) 



determines which of B{T)q^ and D[T)q\ dominates the 
denominators in Eq. (jVII.42p . 
The locus EC , given by 



EC : qz = (C)"' , 



(VII.69) 
(VII.70) 



is where we crossover between the random compression 
and the random field in Eq. (|VII.42p . while AC, given 
by 



AG: qs^iC) , 



(VII.71) 



The above result only applies for both qs^z each much 
larger than its associated inverse correlation length 
Once this ceases to be true, there will, as in the case of a 
generic direction of q just discussed, be two subsequent 



plays the same role between the random tilt and the ran- 
dom field. 

In these six regions, the if-dependence of (|u((7)p) is 
given by 
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2^?/^ — T?t — 6 

ilzS.z) '"^ , region 1 

(gsCf ) , region 2 

[j^) {qzD , region 3 
, region 4 

[j^J [l^iz) > regions 

[iTj [Isis) 1 region 6 



(VII. 72) 



where, just to summarize everything in one place, the This in turn imphcs for the hght scattering: 
boundaries of the six regions are given by: 



FG : 




(VII.73) 


EF : 


Qz = C\ 


(VII.74) 


DF : 




(VII.75) 


OF : 


Qz = -ris^ 

t,z 


(VII.76) 


BC : 


Qz = iCz)-' , 


(VII.77) 


AC : 




(VII.78) 



Qs(t-) (9^6) , region 3 



I{q) cx { 



(€)■ 



(6 



region 1 
region 2 

region 3 
region 4 



-3 



" (iziz) regions 

/j«x»7K + |r)B-2 3 _ . 

i^j (^") 9s ' region 6 



(VII.79) 



The two characteristic lengths and arc obtained 
by equatin g the random field piece D{T)B{T)-i{T)) 
of (|VII.42p to the random tilt and random compression 
pieces, in the appropriate regime of q, and are given by 



with 



Cz 



N 



= 3 - ^7]K + Vt- 



(VII.80) 
(VII.81) 

(VII.82) 
(VII.83) 



Now imagine experiments in which one varies the light 
scattering wavevector q along the two different experi- 
mental loci we have indicated in Fig. U) 



Along locus (1) we fix at and vary q^. We predict 
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His) cc { [irj 



9. 

icr' « qs « c 



(VII.84) 



Thus, a log-log plot of I(qs) verses Qs will be a straight 
line with slope 2riK — r]t — 4 « —3.5 for qs » ^7^, 
where in the approximate equality we have used our e- 
expansion estimates (jl.lOp and (jI.12p for 77i<- and 774 . This 
straight line is then replaced by a horizontal section for 
^ 9.5 ^ Cs~^7 followed by another straight section 
with slope —1 for smaller qs 's. We can read off £,s and 



from the place where the slope changes. Fitting ^ to 
[T — Tac) determines Vs- Knowing and allows 
us to calculate from Eq. (|VII.80|) . 

Along locus 2 wc vary q^ and keep qs fixed at a value 
in the range <C (/s ^ We find 



I{.qz) oc < 



{17) (9-^^) 



9^ > 

-3 



(fj) 9, « q. « (ei)- 



(VII.85) 



,-3 



92 < 



Using an approach similar to that just described for locus 

I, we can obtain {2riK — tit — 4:)/(^z, £,z, CIj ^-nd hence v^, 
and . These results combined with those from locus 1 
allow us to find riB,K,t- Then plugging Vg and r]B,K,t into 
Eqs. (|1.211 II.19|) gives rj^ and f/j, respectively. 

Finally, since rjB.K.t, i^s,z, /?, and a can all be measured 
experimentally, the exact scaling relations Eqs. p.l7| 

II. lSi II.20[) can be examined. 

In the C phase the relation between the light scatter- 
ing intensity and the fluctuations is given by a formula 
similar to (|VII.39P : 

iiq)^Y.^'^AA'i} (vn.86) 

ij 

where 

d, (q) ^ miwW)- (VII.87) 

The coefficients E'^^ are independent of q. However, their 
value is affected by the tilt angle of the layer normals, 
and therefore has a weak dependence on the tempera- 
ture. Since we have argued in Sec. |Vl] that {\u'{q)\'^) 
has the same expression as we conclude that 

the C phase has qualitatively the same light scattering 
behavior as the A phase. 



VIII. STABILITY OF THE SMECTIC PHASE 
AT THE CRITICAL POINT 

All of our discussion so far has assumed that the smec- 
tic state remains stable right up to the AC critical point. 
There are two conceivable mechanisms by which this as- 
sumption might fail: first, the orientational order of the 
smectic layers could be destroyed; and second, the in- 
tegrity of the smectic layers themselves could fail via the 
proliferation of dislocations in the smectic layers. 

In this section, we argue that neither of these hap- 
pens, and that the smectic state is stable right up to, 
and through, the AC transition. 

Consider first the orientational order of the system. In 
deriving our model we have implicitly assumed that the 
both the real space fluctuations {\5n{r)\'^) of the ncmatic 

directors nir), and those (|(5iV(f)p) of the normals N{r) 
to the smectic layers, are small, so that the expansions 
n = z + Sn and N = z + SN are justified. Thus, for 
our theory to be valid, long-ranged orientational order 
must exist. To show that it does, even at the AC critical 
point, at least for sufficiently small disorder, we need to 

calculate {\Sn{r)\^) and {\SN{f)\^). 

We can obtain both of these quantities by calculat- 
ing (|V±itP), since, as shown in section ITTl fluctuations 
of both SN^r) and Sn(f) away from V^w are massive. 



23 



Hence, up to small corrections, (|(5n(r)p) = (|(5A^(r)p) = 
(IV^^P)- 

The most problematical point is right at the AC tran- 
sition, since there the fluctuations are biggest. Hence, if 

(|Jn(r)p) and (|(5Af(r)p) are finite at this point, they'll 
be finite everywhere. 

We will now show that they are, indeed, finite right at 
the AC transition. 

Because of the anisotropic scaling, is the 

dominant part of (|V±Mp). Using our earlier equation 
(|VlL9|l for (|u(g)|2), we obtain: 

1 r Mq)Q^s 



(27r)3 



with 



cx L 



v 



r' 



(VHI.l) 



-I + Tit + riK + 



2 ' 



(VHL2) 



where L is the infrared cutoff in the s directions. Thus, 
r' < is required for the convergence of (|VuP) and 
hence the convergence of ((Jn,;)^). Clearly, in d = 3 
our numerical estimate to first order in e does satisfy 
this condition. This shows that at the AC transition 
point point, the smectic is stable against orientational 
fiuctuations. 

In previous sections we have shown that right at the 
critical point the long-ranged translational order is de- 
stroyed in the presence of the quenched disorder. Thus 
one might reasonably question what is the difference be- 
tween the smectic and the isotropic (nematic) phases, 
and even worry whether the smectic phase survives at 
the AC critical point at all. However, there are examples 
in which the low-temperature and the high-temperature 
phases are not distinguished by the order of the system, 
but rather, by whether the topological defects are bound 
or unbound. The most famous one is the Kosterlitz- 
Thouless transition in the d = 2 XY model [l^l • Another 
famous example is the " Brag g glass" in pinned supercon- 
ducting fiux lines [3, [3ll - l33| . We will now show that the 
same thing happens in our problem. Since similar calcu- 
lations have been given in Refs [H, in great detail, 
and ours are virtually the same, our description will be 
very brief. 

The starting point of the theory is the tilt only har- 
monic Hamiltonian at the critical point: 



H 



d-^r 

7 



K 



2 



B 



h ■ Vw 



(vni.3) 



where the random field disorder is not included. It can 
be justified that the random field disorder is irrelevant 
when dislocations are included. 



When the smectic contains dislocations, the displace- 
ment field u is no longer single valued. Mathematically 
this can be represented as 



V X VM : 



(VIII.4) 



with 



aM,t{s,)d^ (f - f,{s,)) ds, , (Vin.5) 



where Si parameterizes the path of the z'th dislocation 
loop. Mi is an integer giving the charge of that loop, 
t{si) is the local tangent of the loop, and is a position 
of the point on the loop parameterized by Si the loop. 
Furthermore, equation (|VIII.4p implies 



V • il/ = 0, 



(vin.6) 



which means that dislocation lines can not end in the 
bulk of the sample. 

To obtain a dislocation Hamiltonian, we need to trace 
over field u which is constrained by equation (jVIII.4p . 
This can be done in the following standard way. We 
separate the field v = Vu into 



V = Vd + 5v, 



(vin.7) 



where Vp minimizes Hamiltonian (|VIII.3|) for a given dis- 
location configuration TO(r), 5v can be viewed as the fluc- 
tuation from the ground state. Inserting (|VIII.7I) into 
Hamiltonian (|VIII.3|) . we find that Vp and 5v are decou- 
pled due to the construction that Vd minimizes Hamil- 
tonian (|VIII.3|) . Thus we obtain the effective model for 
fri{r): 



Hd 



E 



2BT, 



{K^ql\mM? + BKql\mh{q)[' 



+B^\ms{q)\'') + rn{q) ■ d{q)\ , 
where Tg = ql + / B)ql + Xlq^, and 



(VIII.8 



qxh 



+ 



jqhh + Kq^sj x q 



q- h 



BTgq' 



(vin.9) 



By putting the model on a a simple cubic lattice, the 
partition function can be written as 



-S[m] 



where 



TO(r) 
S[m] 



— K(f),nj,(f),n^(f)], 



1 

T 



(VIH.IO) 



(VIH.ll) 



(VIII.12) 
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where the n^'s are integers, d is the cubic lattice constants 
used in the discretization, and |m(r)p is the core 

energy term coming from the core of the defect hue that 
is not accurately treated by the continuum elastic theory. 

To cope with the constraint V • m = 0, we introduce 
an auxiliary field 9(f), rewriting the partition function 
equation (jVIILlOp as 

Z^llJ d9{r) e-^["l+*^^^''«^■™«'''/^(VIII.13) 

r m(r) 

where the sum over m(r) is now unconstrained. 

Then we introduce a dummy vector field A to mediate 
the long-range interaction between defect loops in the 
Hamiltonian (|VIII.8[) . This is accomplished by rewriting 
the partition function as 



^ = n / dOif)dA{f) J2 e 

r m(r) 



-S\m,e,a] 



SiV-A)5{A,) 



(VIII.14) 



with 



-i—S/e{r) + d^[iA{r) + d{r)] 



-E 



( 



2T ^ \ K-/q] 

q ^ 



1 



(VIII.15) 



The sum over m(r) is just the "periodic Gaussian" 
introduced by Villain (33 |. Then the partition function 
equation (jVIII.14[) can be rewritten as 



Z = YlJ d9{r)dA{f) (5(V • A)S{A,) 



X exp 



(f + Xi) 



^lA,{q)\^ + ^\AM\' 



Kg] 



7 



(VIII.16) 



where the 27r-period Villain potential Vp{x) is defined as 
e-^p(-) = J2 6-"'^=/^+"". (VIII.17) 

— oo 

Since Vp{x) has sharper minima for smaller Ec/T , which 
corresponds to higher temperature, raising the tempera- 
ture in the original model is equivalent to lowering tem- 
perature in the dual model equation (jVIII.16[) . 

Standard universality class arguments imply that the 
model equation (|VIII.16p has the same universality class 



as the "soft spin" , or Landau-Ginsburg- Wilson, with the 
complex "action" 



r.a. '- ^ ■'^ ^ 

d) + + udiVoii + E 



Kq. 



(VIII.18) 



where ip is a complex order parameter whose phase is 
0{r). Because of the duality transformation's inversion 
of the temperature axis, the reduced temperature td is a 
monotonically decreasing function of the temperature T 
(of the original dislocation loop model), which vanishes 
at the mean-field transition temperature of the disloca- 
tion model (jVIII.16[) . Disorder is included in the model 
(|VIII.18[) through a(r), which is related to the random 
tilt field h{r) by equation (jVIILOp . 

Because of the duality inversion of the temperature 
axis, the ordered phase of the dual model (|VIII.18p cor- 
responds to the disordered phase (dislocation loops un- 
bound) of the original dislocation model. 

An complete analysis of the dislocation loop unbinding 
transition described by the model (jVIII.18|) is beyond 
the scope of this paper. The goal here is to know if a 
dislocation bound state ever exists near the critical point 
of the AC transition. Let us check the one-loop graphical 
correction to the dual temperature: 



tR — to 



caH^ f d^q ^t{B^ql+l^ql) 



T2 



(27r)3 



.(VIII. 19) 



If K and Aj are treated as constant it is easy to show 
that the integral diverges. This result, if true, would 
imply that the smectic phase becomes unstable against 
dislocation unbinding sufficiently close to the putative 
AC transition. 

However, this conclusion only holds within the har- 
monic approximation. In sections IIIII and IIVI we have 
shown that anharmonic effects are important. A crude 
way to include their effects is to replace the constant K 
and A in the harmonic calculation above with the anoma- 
lous, wavcvector-depcndcnt K{q) and At(q) found earlier 
in section IIVI Implementing this in the integral and re- 
quiring it to converge, we obtain a restriction on r//^ and 



1 



(VIII.20) 



This condition is satisfied by the values oi rjs and r]t 
obtained by the e-expansion to 0(e) in the physical di- 
mension d ~ 3. This implies that the the second-order 
smectic A-C phase transition is stable against dislocation 
unbinding. 
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IX. CONCLUSION 



We have developed a theory of the smectic-A to 
smectic-C phase transition in a biaxial disordered envi- 
ronment. Using an e = I — d expansion, we have shown 
that the phase transition can be second order. The crit- 
ical exponents were calculated to first order in e. We 
also calculated the smectic correlation functions near the 
critical point, which lead to exotic X-ray and light scat- 
tering patterns. All these predictions should be testable 
provided that the biaxial disordered environment can be 
realized in experiments. 
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Appendix A: 



There are two ways of continuing our model Eq. (jll.ll[) 
from three dimensions to higher spatial dimensions. In 
the main text we have focused on the "hard" continu- 
ation, in which we keep the /i-space d — 2 dimensional 
and the s-space one dimensional. In this appendix we 
briefly describe the calculation which uses the soft con- 
tinuation, in which we keep the /i-space one dimensional 
and the s-spacc d — 2 dimensional. 

The RG procedure is essentially the same as that used 
in the hard continuation. After integrating out the free- 
dom of the fast varying part of u{r) and rescaling, we 



obtain the following RG flow equations: 

[d-ujh + t^z + 2x- 2] 7, 



c?7 
H 
dB 

u 

dK 

m 

d£ 

d£ 

dw 



dAt 

de 

dD 



+ cj/i - + 2x - 2 - igi 
d + + uj^ + 2x - 6 + -gi 



d + ujh + 3x- 4:+ igi - 



3 

rf + a;;, + W3+4x-6--^ 
8.92 



5. 9. ■ 



w, 



d + u!h+uiz + 2x- 4:+ —gi 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



(A.5) 
At, (A.6) 



d + Uh + ^z + 2x-4 + -{gi-h) D 
-lK{g2~~gil (A.7) 



where gi and 52 are two dimcnsionless couplings defined 
by: 



51 = Ca-2{9lBfAt{BhY'^K-^k''-^ 

52 = Cd-2{w/B)At{Bhf'^K-^K''-^ . 

These couplings gi and 52 flow according to 

. 2^ 



dgi 
de 
dg2 



£+3^1 -32, 

e- a 0.92 + ^31 

8 52 8 6 



(A.8) 
(A.9) 

(A.IO) 
(A.ll) 



where e = 4 — d. 

Once again, it is convenient but not necessary to make 
a special choice of Xj w^,, and oJz such that B, K, and 7 
are fixed at their bare values. This requirement can only 
be fulfilled by 



(A.12) 
(A.13) 
(A.14) 







1 




= 2 


3 . 




= 2 
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This choice turns fiow Eq. (jA.7p into 
dD 



g-9i + 8 (-91 -92) 



D + -K{g2 - .gi). 



(A.15) 



Flow Eqs. (|A.7IA.10IA.11[) have three fixed points in 
total. One of them is the trivial Gaussian fixed point 
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gl = g2 = 0, D = 0. Another one was also found pre- 
viously for the 7Ti = 1 Bragg Glass probleni[l3, which 
preserves the symmetry that the Hamiltonian is invari- 
ant under any rigid rotation of the liquid crystal about 
the h axis. This fixed point is: 



9i 
D* 



92 

0, 



(A.f6) 
(A.17) 



Both of these fixed points are unstable in two or more 
directions. 

The third fixed point is: 



9i 

92 

D* 



(A.18) 
(A.19) 
(A.20) 



which is stable in two directions, and unstable in a third. 
Therefore, this fixed point controls the critical behavior 
of the phase transition. 

With this fixed point in hand, the anomalous elasticity 
and the critical exponents can be calculated by following 
the steps described in Sees. llVl andlVl Here we only give 
the results: 
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2 



i'h = f 



= f 



(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 
(A.26) 
(A.27) 
(A.28) 



Now we compare the numerical results obtained from 
the two continuations. In d = 3, to first order in e and e 
we find, for the hard continuation discussed earlier. 



VB 


= 0.2222, 


(A.29) 


Vk 


= 0.2963, 


(A.30) 
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(A.31) 




= 0.6296, 


(A.32) 
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(A.33) 




= 1.1296, 


(A.34) 
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= 0.3889, 


(A.35) 
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= 0.1111. 


(A.36) 




FIG. 7: Graphical correction to Ac 



while for the soft continuation treated in this appendix, 



'TIB 

Vt 
a 



0.1071, 

0.05357, 

0.01786, 

0.5938, 

1.1741, 

1.1473, 

0.2188, 

0.1607, 



(A.37) 
(A.38) 
(A.39) 
(A.40) 
(A.41) 
(A.42) 
(A.43) 
(A.44) 



Unfortunately, the two sets of results are not in good 
agreement. We trust the hard continuation more, since e 
is smaller (e = i) for that configuration than the expan- 
sion parameter e ~ \ for the soft continuation. 



Appendix B: Wavevector-dependences of Irrelevant 
Disorder Variances 



In our RG calculations in section IIIIl we did not in- 
clude in the Hamiltonian (jIII.I[) the random compression 
disorder 



dzUp , 



(B.l) 



because simple power counting shows that this term is 
irrelevant. However, while it has no effect on the crit- 
ical exponents, the disorder invariance Ac itself devel- 
ops strong power-law dependence on wavevector at long 
wavelengths. We need to know this power-law exponent 
to calculate the fluctuations caused by the random com- 
pression, which have a significant effect on the light scat- 
tering intensity. 

The most important one-loop graphical correction to 
Ac comes from the Feynman diagram in Fig. [71 from 



27 



which we obtain 
r.2 r> 



6A, 



5_ 
2 



52 A? r dp. 



2 7-00 27r 
11 At^gi 
512V2 A'A2 



^ d^~^p^ f> dp, 8^ ,^^4 



d-2 



(27r) 



(B.2) 



Thus, using the length and field rescalings we used in 
section Hill the RG flow equation of Ac to one- loop order 
is given by 



di 



2x+{d- 2) ujh-^z + l + 



11 



512V2 



33 



A, 



where the dimensionless coupling is defined as 

BAt 



93 



KArA^ 



51 



(B.3) 



(B.4) 



Combining this RG fiow equation with the flow Eqs. 

(Hiltti Hirn innii ulth), we get 



dgs / 8 11 

-di = ['-r-^''^''' 

which fiows to a stable nonzero fixed point 
V2 



93 



99 



(9216 - 4096e) . 



(B.5) 



(B.6) 



Now the wavevector dependence of Ac can be calcu- 
lated by using trajectory integral matching on the cor- 
relation function {u{—q)){u{q)) ~ [Atql + Acg^)G(g)^, 
Performing a calculation which is essentially the same as 
those in section ITVl we find 

AM = Ac(efg.)""Vc(X,>^) , (B.7) 



with the anomalous exponent 



(B., 



Here fc{X, Y) is another crossover function of the scaling 
variables X and Y given by (|IV.10|) and (jIV.llj) . respec- 
tively. Note that 77c is nonzero even to zeroth order in e, 
which is quite common for irrelevant variables. 



We can also obtain Eq. (jB.Sp through an exact scal- 
ing relation between 77c and other anomalous exponents. 
This scaling relation is implied by the fact that 93 flows 
to a non-zero stable fixed point. For large enough £, 
reaches a fixed point and thus 



d In gz 
d£ 



^0. (B.9) 

This equation, after rewriting 53 in terms of other pa- 
rameters using its definition (|B.4p . leads to 



(B.IO) 



d\nB dlnAf din 51 d\nK din Ac 



di 



de 



de 



Then we replace the graphical part in flow Eqs. pil.lll 
nil. 121 lni.l5IB.3|) by —775, riK,t,ci respectively, and plug 
them and Eq. piI.19p into the above equation. The 
rescaling factors x and uj vanish, and we are left with 



m - Vb - Vk- 



(B.ll) 



Using the already known values of 774, 77s, and rjx to 0(e), 
using ()B.11|) also leads to the epsilon expansion for 77c; 
reassuringly, this expansion is exactly (jB.8[) . 
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